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A continuous sequence of infinitesimal unitary transformations, combined with an operator product ex- 
pansion for vertex operators, is used to diagonalize the quantum sine-Gordon model for (3 2 £ (2ir, oo). The 
leading order of this approximation already gives very accurate results for the single-particle gap in the 
strong-coupling phase. This approach can be understood as an extension of perturbative scaling theory 
since it links weak to strong-coupling behavior in a systematic expansion. The method should also be 
useful for other strong-coupling problems that can be formulated in terms of vertex operators. 
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Perturbative scaling arguments play an important role 
for analyzing a large variety of physical systems with 
many degrees of freedom. For strong-coupling prob- 
lems, however, the perturbative renormalization group 
(RG) equations lead to divergences in the running cou- 
pling constants and the perturbative RG-approach be- 
comes invalid. In condensed matter theory the well- 
known paradigm for this kind of behavior is the Kondo 
model: The perturbative scaling equations still allow one 
to identify the low-energy scale of the Kondo model, but 
by themselves they do not lead to an understanding of 
the physical behavior associated with this energy scale. 
Wilson's numerical RG Q could remedy this problem, 
but an analytical RG-like approach that links weak to 
strong-coupling behavior in an expansion that can be 
systematically improved would still be desirable for many 
strong-coupling problems. 

In this Letter it will be shown exemplary how Wegner's 
flow equations Q can provide such an analytical descrip- 
tion for a weak to strong-coupling behavior crossover. In 
the flow equation approach a continuous sequence of in- 
finitesimal unitary transformations is employed to make 
a Hamiltonian successively more diagonal. Large energy 
differences are decoupled before smaller energy differ- 
ences, which makes the method similar to the conven- 
tional RG approach. However, degrees of freedom are 
not integrated out as in the RG but instead diagonal- 
ized. A similar framework that contains Wegner's flow 
equations as a special case has independently been devel- 
oped by Glazek and Wilson (similarity renormalization 
scheme) ||. 

The model under investigation in this Letter is the 1 + 
Id quantum sine-Gordon model pj 



ut 2 cos \J3(f>(x)\ 



with the commutator [II(a;), 4>(x')] — —i6(x 



Regu- 



larization with a UV-momentum cutoff A cx r is im- 
plied and u, (3 > are dimensionless parameters. The 



sine-Gordon model is one of the best studied integrable 
models and it has been solved using the inverse scatter- 
ing method ||. This model is therefore a good test case 
for studying the new approach. We will be interested in 
the universal low-energy properties (E <C A) for small 
coupling constants u. It should be emphasized that the 
integrable structure underlying the inverse scattering so- 
lution will not be used in the approximate flow equation 
solution; the new method can also be used when non- 
intcgrable perturbations are added. 

The phase diagram of the sine-Gordon model consists 
of a gapped phase for 1 < 8ir with massive soliton ex- 
citations and a gapless phase for 1 > 87r with mass- 
less solitons Q . The phase transition between these two 
phases for f] 2 /8n = l + 0(u) is of the Kosterlitz-Thouless 
type. In the massive phase the perturbative scaling equa- 
tions H lead to an unphysical strong-coupling divergence 
of the running coupling constant u. The inverse scatter- 
ing solution m furthermore shows the existence of bound 
soliton states (breathers) in the spectrum for /3 2 < 4ir 
while such bound states are absent for (3 2 > 4ir. For 
/3 2 = 47T the sine-Gordon model can be mapped to a 
noninteracting massive Thirring model , which in turn 
can be diagonalized easily leading to the identification of 
the massive solitons with the Thirring fermions The 
sine-Gordon model is also related to a variety of other 
models like the spin-1/2 X-Y-Z chain, a Id Fermi sys- 
tem with backward scattering and the 2d Coulomb gas 
with temperature T — f3~ 2 and fugacity z oc it, the IR- 
unstable fixed point corresponds to T — 1 /8ir || . 

It will be shown that the flow equation approach gen- 
erates a diagonalization of the sine-Gordon Hamiltonian 
both in the weak-coupling and in the strong-coupling 
phase in a systematic expansion that can be successively 
improved: No divergences of the running couplings are 
encountered in the strong-coupling regime for 1 > 2ir. 
The soliton mass is found to be in very good agreement 
with the inverse scattering solution. The crossover from 
weak- to strong-coupling behavior can be described and 
using this the soliton dispersion relation for example can 
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be analyzed on all energy scales. 

For the purposes of this Letter it will be more con- 
venient not to use a regularization of the sine-Gordon 
model with an explicit momentum cutoff, but instead to 
"smear out" the interaction term 



H = dx 



I 



Ii(xf 



(1) 



P / dy c(y) <f>(x + y) 



with the Lorentzian c(y) — a/(2ny 2 + na 2 /2) S(y) 
fiof . This does not affect the universal properties for 
small energies B < A a a' 1 . We expand the fields in 
normal modes 



4>(x) 



n(aO 



4 -S 



-ikx 



(<Tl(k) + <T 2 (k)) 



\k\e~^ x (ax(fc) - a 2 (k)) 



where J2k ^ 17 J2n=-oo w ^ n & = 2irn/L. L is the 
system size. The basic commutators are (k, k' > 0) 
[cri(-fc),cTi(/c')] = [a 2 (k),a 2 (-k')] = 5 kk , L/2tt, notice 
aj(—k) = a i (k). The vacuum is defined by a%(— k)\Ct) = 
a 2 (k)\Q) = for k > 0. 

The flow equation approach |g| (see also jll)) gener- 
ates a family of unitarily equivalent Hamiltonians H(B) 
as a function of a flow parameter B (with dimension 
(Energy) -2 ), where H(B = 0) is the initial Hamiltonian 
and H(B = oo) the final diagonal Hamiltonian. This 
flow is generated by the differential equation 



dH{B) 
dB 



[r)(B),H(B)] 



(2) 



with rj(B) = —rtiBy some antihermitean generator. 
Generically, Eq. (||) leads to the generation of new in- 
teraction terms not contained in the initial Hamiltonian. 
We therefore write H(B) = H + H int {B) + H ncw (B) with 



H 



p>0 



t p (aiO)o-i(-p) + o- 2 (-p)a 2 (p)) 
) 

Hi n t(B) = J dxi dx 2 u(B; x\ - x 2 ) (3) 

x(Vi(a;xi)V 2 (-a;x 2 ) +h.c.) . 

Here Vj(a;x) are normal ordered vertex operators with 
scaling dimension a{B) d = /3(B) /y/Zn 

s., 

, p ' 



Vj(a;x) 



exp f ± 



-f \p\-ip- 



where + (upper sign) for j = 1 and — (lower sign) 
for j = 2. To avoid confusion the initial values of 



the couplings will from now on be denoted by uq,Pq. 
H(B = 0) is identical to the sine-Gordon Hamilto- 
nian (@) for H ncw {B = 0) = 0, a(B = 0) = Aj/VJtt 
and u{B = 0;x) = u 6(x) {2na/ L) Q ( B =°) 2 /2vra 2 . 

Wegner's idea for constructing a suitable generator ij 



is to choose r](B) =' [H ,H int (B)] gj. This gives 

7] = -2i J dxdy (d y u(y))(Vx{a; x)V 2 (-a; x - y) + h.c). 

Using this generator, matrix elements connecting states 
with large energy differences are eliminated for small B 1 
while matrix elements coupling more degenerate states 
are eliminated later for larger values of B. 

First we evaluate [r],Ho]. This leads to du/dB = 
4 d 2 u/dx 2 , which makes the interaction increasingly non- 
local along the flow due to the decoupling procedure. In 
Fourier components u(B;x) = ^2 p Up(B)e~ lpx one finds 



u p {B) 



u{B) 
4:TT 2 a 2 



-4p 2 B 



2ixa 



u{B) will turn out to be the running coupling constant of 
the flow equation approach and remains finite also in the 
strong-coupling phase. u{B) is dimensionless, initially 
u(B = 0) = uo- The term [77, Hint] leads to new interac- 
tions that have to be truncated to obtain a closed set of 
equations. The approximation used here is to take only 
operators with small scaling dimensions into account, i.e. 
to neglect more irrelevant operators. It is generated by 
truncating the operator product expansion (OPE) of two 
vertex operators in the following way 



Vj{a;x)Vj(-a;y) = — 



1 



2V [a T i{x-y)r 
I T ia(x - y)Y,V\p\e-> lpl - ipx *j(p) + ■ 



(4) 



p#0 

The approximation can be systematically improved by 
going to higher orders in this OPE. The term [77, Hint] 
contains commutators with the structure 

[Vt(a; Zl )V 2 {-a; z 2 ), V 2 (a; z'M-a; z[)\ (5) 
= -{^(o; z 1 ),V 1 (-a; z[)}V 2 (a; z' 2 )V 2 (-a; z 2 ) 
+Vi(a; zi)Vi(-a; z[){V 2 {a; z' 2 ), V 2 (-a; z 2 )} 

and terms where a — > —a in one argument of the com- 
mutator (^) . After normal ordering, the latter terms lead 
to interactions Vi(2a; Zi)V 2 (— 2a; z 2 ) with larger scaling 
dimensions. These will be neglected, but the terms gen- 
erated by (^) will be included. 

For a = 1 (f3 2 = An) the vertex operators de- 
scribe fermions and the OPE (Q) to all orders gives 

{Vj(—l;x),Vj(l;y)} L5(x — y). Since this is a c- 
number, no higher order interactions are generated in (0) 
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and the flow equations close. The flow equations there- 
fore recover the equivalence of a sine-Gordon model with 
01 = 4tt to a massive noninteracting Thirring model |t]] 
and readily diagonalize the latter. 

In general we evaluate (||) using (Q). The dominat- 
ing contributions decaying most slowly with B can be 
identified in closed form [Q. Two structurally different 
interaction terms are generated: One term contributes 
to ffnew and is discussed below (see Eq. (|l0|)). The 
other term has the structure a\(k)(j2(—k). Integrating it 
from B to B + dB, one generates J2 k Wk\k\ai(k)a 2 (— k) 
with infinitesimal coefficients Wk- This new interac- 
tion can be removed by a further infinitesimal unitary 
transformation with the structure e~ u H(B)e u , where 
U = ^X)fc>o u,fc ( <Tl (^) cr2 (~ k) — h.c). This new infinites- 
imal unitary transformation yields a modification of the 
scaling dimension of the vertex operators in H mt and a 
flow of the coupling constant u(B). In terms of a loga- 
rithmic dimensionless flow parameter £ = I In (32Ba~ 2 ) 
one derives ll2l 



d{3 2 {£) 
dl 



47rr( - 1 + (3 2 (l)/An) 



(6) 



x exp Al 



dt' p 2 (e') 



The running coupling constant flows according to ii(£) 
uq exp (F(£)j, where 



F(£) 



£(3 2 (l)- / d£'[3 2 (l') 



(J) 



In the strong-coupling phase the flow terminates at 
(3 2 (oo) — An due to the divergent T-function in 
(^). In our approach 1 = An is therefore an at- 
tractive strong-coupling fixed point. For comparison 

with the RG-equations || one can introduce u(£) 

and rewrite (p[) as two cou 



dcf 



uq exp 
pled differential equations 

d{3~ 2 (£) 



dt 

du(£) 
dl 



AnT(- 1 
P 2 (£ 



= 2- 



An 



p 2 (£)/An 

u(£) 



■u 2 (£) 



(8) 



with P(£ = 0) = Po,u(£ = 0) = uq. For (3 2 = 8n 
Eqs. (||) coincide with the two loop scaling equations ||: 
Depending on the value of itn, the sine-Gordon model for 
[3 2 > 8n flows to either f3 2 (00) — An (strong-coupling) or 
P 2 (oo) > 8n (weak-coupling). Eqs. (^) therefore repro- 
duce the Kosterlitz-Thouless phase diagram. Also the 
hidden SU(2)-symmetry in (0) for /3g = 87r(l±u ), u < 
1 is recovered although our approximation scheme does 
not manifestly respect this symmetry ]l3| ]. 
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FIG. 1. Soliton mass as a function of the coupling constant 
for various values of The full lines are constrained fits 
of the power law behavior am oc u^ 2 /3 °' /4 ' r ^ to the flow 
equation results (open circles) with the proportionality con- 
stant being fitted. The dashed line is the case 0q = An where 
the flow equation approach agrees trivially (see text). 

In the strong-coupling phase a gap opens in the spec- 
trum and the low-energy excitations are fermionic || . In 
our approach this follows most easily by approximating 
(|) for large B (such that \P 2 (B) - An\ < 1) as 



H int (B) 



dx\ dx: 



u(B) 



■. exp 



(^1(^1)^2(^2) 



-^2(^2)^1(^1 



(#i - x 2 f 



165 



(9) 



dcf 



with fcrmions ipj( x ) = Vj'( — l) x )- The running coupling 
constant u(B) approaches a finite value in this limit |l4) . 
The asymptotic value u(oo) can therefore be interpreted 
as the renormalized coupling constant parametrizing a 
quadratic Hamiltonian Hq + H lnt (B) that describes the 
low-energy behavior of the initial sine-Gordon model. 
The RG strong-coupling divergence of the running cou- 
pling constant is avoided since the interaction ([)]) in this 
effective low-energy Hamiltonian becomes increasingly 
non-local for B — > 00. The gap A = 2m in the spec- 
trum is obtained in terms of the renormalized coupling 
constant u(oo) that sets the mass of the effective low- 
energy Thirring fermions with m — u(oo)/a. 

The soliton mass derived from Eqs. (^||) can be 
compared with the predictions of perturbative scaling: 
E.g. for [3 2 — 8tt(1 — itn) (corresponding to the hid- 
den SU(2)-symmetry) one can show analytically m oc 
u5exp(-l/2w )/a for u J. with r = (l + 7)/3 w 0.526 
(7 w 0.577 is Euler's constant). This agrees well with 
the third order RG result trq = 1/2 and the same expo- 
nential dependence For smaller values of /3n one 

finds the power law behavior m cx u^/j 2 /3 °^ 4 ^ /a known 
from the inverse scattering solution B : This is very well 
confirmed by the numerical integration of the flow equa- 
tions (Fig. 1). The proportionality constant depends only 
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on fig for uq I 0, in particular m — uq/o, for 1 = Air as 
known from the noninteracting Thirring model [Q . 

The flow equation solution not only provides the exci- 
tation gap, but one can also obtain information about the 
crossover, e.g. the full dispersion relation: The Hamilto- 
nian for B = oo takes the form H(oo) = Hq + H nev! (oo) 
since the interaction term Hi nt (B) is eliminated. One 
can verify that H new (B) does not modify the flow of 
(3(B) and u(B) as derived above The new terms 

generated during the flow in H nevr (oo) can be split up as 
ffncw(oo) = i/dia g (°o) + H Tes (oo), where H dias (oo) con- 
tains the terms that follow in leading order of the OPE 
from [77, ^Hint] in (||), while H ICS (oo) formally contains ev- 
erything not taken into account in the present order of 
the OPE. Integration of the flow equations gives Jl2| 

ffdiagM = J2^)(pHp)Pi (p) + Pi (-p)pH-p) 

p>0 

+pU-p)p 2 H) + P2 (p)Pkp)) (io) 

with certain coefficients u) p (oo). Here Pj{p) = 

J dx e- ipx V 3 (-a(B p );x) and B p d = l/Ap 2 . The spec- 
trum can be analyzed easily since [Hq, i?diag(oo)] = 0: 
In leading order the single-particle (soliton) excitations 
of H + ffdiagM are Pl(k)\Cl) for k > and P 2 f (fc)|0) 
for k < 0. The single-hole (antisoliton) excitations are 
Pi(k)\n) for k < and P 2 (k)\n) for k > 0. In the strong- 
coupling phase for /3q > 4-7T the resp. excitation ener- 
gies are very accurately (but not exactly) described by 
El = k 2 + (u(oo)/a) 2 in the small coupling limit: There 
are /^-dependent universal corrections in the crossover 
region k = 0(m) that vanish for f3 2 — > An and reach at 
most of order 2% for 0$ — 8tt. The character of the ex- 
citations varies from scaling dimension a(B = 0) to the 
low-energy Thirring fermions with a(B = 00) = 1. In 
the weak-coupling phase the spectrum remains gapless 
and Eh — \k\ for k — * 0. Notice that the elementary 
excitations are expressed with respect to a transformed 
basis since H(oo) and H(0) are related by a complicated 
unitary transformation. 

For 0Q < 2ir the differential equations for lo p (B) lead 
to divergences since then the cos(/?o'/ , (^))-perturbation 
is too relevant, limiting the approach to 0$ > 2tt. For 
2ir < 0q < Air the single-particle spectrum is still well 
described by (^0|), our approximations become better as 
00 t 47r. Higher orders in the OPE are nevertheless re- 
quired to study the formation of bound states for /3q < Att 
H due to residual interactions in H ICS (oo). 

Summing up, we have applied a continuous sequence of 
infinitesimal unitary transformations, combined with an 
operator product expansion for vertex operators, to the 
quantum sine-Gordon model with 1 £ (2tt,oo). The 
approximations are systematic since more terms in the 
OPE can successively be taken into account and will 
not endanger the stability of the strong-coupling fixed 



point. The results for the soliton mass in the strong- 
coupling phase agree with two loop scaling predictions 
for 1 w 8ir (approximate agreement was even found to 
three loop order) and exact methods applicable for 
smaller 1 . The full dispersion relation could be obtained 
and the crossover from weak to strong-coupling behavior 
described. The method also allows one to study corre- 
lation functions |l6| ] and other strong-coupling problems 
that can be formulated in terms of vertex operators. 
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